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ABSTRACT 

We have computed seismic images of magnetic activity on the far surface of the Sun by using a 
seismic- holography technique. As in previous works, the method is based on the comparison of waves 
going in and out of a particular point in the Sun but we have computed here the Green's functions from 
a spherical polar expansion of the adiabatic wave equations in the Cowling approximation instead of 
using the ray-path approximation previously used in the far-side holography. A comparison between 
the results obtained using the ray theory and the spherical polar expansion is shown. We use the 
gravito-acoustic wave equation in the local plane-parallel limit in both cases and for the latter we 
take the asymptotic approximation for the radial dependencies of the Green's function. As a result, 
improved images of the far-side can be obtained from the polar-expansion approximation, especially 
when combining the Green's functions corresponding to two and three skips. We also show that the 
phase corrections in the Green's functions due to the incorrect modeling of the uppermost layers of 
the Sun can be estimated from the eigenfrequencies of the normal modes of oscillation. 
Subject headings: Local helioseismology, magnetic activity 



1. INTRODUCTION 

The idea that the observed wavefield at the solar sur- 
face can be used to estimate the wavefield at any other 
location in the Sun has resulted in the development 
of several local helioseismology meth ods such as th e 
seismic- holography technique (Lindsey & Braun 1997). 
Although the technique has bee n used to ch aracterize 
fiows in the solar subphotosphere (|Braun et al.ll2 QQ4) and 
to investigate seismic sources from fiares (|Donea et al.i 
[2OO6), it s application to map areas of strong magnetic 
field (^Lindsev fc Braunlll9 98l has remained a major fo- 
cus. In particular, the mapping of act ive regions in the 
non- visible (or far side) of the Sun (iLindsev & BraunI 



ll99Ql . l2QQQal : lGonzalez Hernandez et al.ll2QQ9f ) has proven 
to be a powerful tool for space weather forecasting. A 
general review on local helioseismology can be found in 
iGizon &: Bird] (j2QQ5). 

The use of phase-sensitive seismic holography to map 
areas of strong magnetic field is based on the fact that 
there is a travel delay or phase shift between waves en- 
tering and e xiting an active region (|Braun et al.l Il992[ 
iDuvall et al.l il996). The technique compares theoreti- 
cal waves going in and out of a particular point in the 
Sun (the focus) with the observed wavefield, ^ ohs(S^-,'^)i 
in some portion of the solar surface (the pupil 7^±), by 
computing the ingression and egression functions given 
by, 

(1) 
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where tq is the location of the focus, R the solar radius, 
Q = (^,^), g -\-{ro, to] R^Q^t) is a Green's function that 
represents the theoretical disturbance at (i?, ^1, t) result- 
ing from a unit impulse origin ating at (vo.to) and Q 
is the time-reverse counterpart fi Lindsev Braunll2QQQbl 
^200^. 

In the quiet Sun, the wave field can be well represented 
by a superposition of gravito-acoustic Green's functions 
of a standard reference model and the ingression and 
egression functions are very similar. However, when the 
focus is located in an area of high magnetic activity, the 
measured signal deviates from that, effectively introduc- 
ing a phase shift between the waves going in and out. 
This can be detected by computing the correlation that 
in the Fourier domain becomes the product of the Fourier 
transforms, i^+(ro,a;) and i^^(ro,cj). 



C(ro) 



if+(ro,w)iJ* (ro,w)dt<j . 



(2) 



The phase of the correlation, ^(ro) = argC(ro), is re- 
lated to the perturbed travel time by (/) :^ ujAt. The 
seismic- holography maps are representations of ^(ro) as 
a measurement of the disturbance of the wave field due 
to the presence of magnetic regions. 

Up to the present, far-side seismic holography has used 
Green's functions calculated using the (acoustic) ray the- 
ory, e.g. [Lindsey & Braun (2000b). Formally, this ap- 
proximation is only valid for high frequencies and high 
angular degrees. In this paper, we develop the formalism 
to determine the Green's functions through their polar 
expansion. Although the method allows a full numerical 
computation both in the radial and angular variables, in 
the present work we have used the asymptotic approxi- 
mation for the radial dependencies of the Green's func- 
tion. On the other hand, the angular factors are com- 
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puted without any asymptotic approximation and hence 
the Green's functions are expected to be more accurate 
for low and intermediate degrees. 

The outhne of the paper is as fohow. In section 
2 we describe the basic wave equations. Section 3 
presents the calculation of the Green's functions using 
the gravito-acoustic ray-path approximation and section 
4 the Green's functions using the spherical polar expan- 
sion. The phase shift introduced by the non-adiabatic 
sub-surface layers of the Sun is reviewed in section 5. 
In section 6, we present an example of the improvement 
achieved by using the new Green's functions to calculate 
far-side maps of solar activity. 

2. BASIC EQUATIONS 

We start the analysis with the adiabatic wave equation 
in the Cowling and local plane-parallel approximations. 



In terms of the scalar field ^(t,r) 



-P 



-1/2 



where p 



is the density of the background state and 5p{t^ r) the 
Lag!:ran^i an pressure fl uctuation, this equation is given 
by — see lGoughl (|1993[ )— 



N^Vl^ = (3) 



where c and N are the adiabatic sound speed and buoy- 
ancy frequency, deflned as usual, represent the terms 
involving horizontal derivatives in the operator and 
the critical acoustic frequency is given by 



^(l-2r.Vff) 



1/2 



(4) 



Here H is the density-scale height and r is an upward 
directed unit vector. 

If the background state does not depend on time, as we 
will assume, the scalar fleld can be Fourier decomposed 
such that the wave fleld is a superposition of monochro- 
matic waves given by ^(r, t) = ^(r, cj) exp (— zco't) and 
Eq. (|3j) becomes 
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We deflne the operator 



Co 



A/'2 1 



cO) 5 



(5) 



(6) 



where A^o , cq and ujco are the buoyancy frequency, sound 
speed and cut-off frequency of a spherically-symmetric 
reference model. In this work we w ill use Model S from 
IChristensen-Dalsgaard et al.l (| 19961 ) . We assume that the 
actual solar oscillations satisfy a wave equation of the 
form 

jC^ = jCo'^^jC'^ = 0, (7) 

with >C'^ accounting for all departures from Eq. ^ (e.g. 
non adiabatic effects) and from the reference model (e.g. 
non spherical perturbations to the sound speed). 

2.1. Integral representation 

Let us take a target point tq, the focus, anywhere in the 
Sun. We introduce Green's functions, G(r|ro;cj), that 
satisfy the source-point wave equation for the standard 
solar model, 



CoG{r\ro;uj) = -S{r - tq) . 



(8) 



Multiplying Eq. ^ by G, Eq.© by ^, subtracting both, 
integrating over a volume V that containing the focus tq 
and, flnally, using the Green theorem, the above equation 
can be written as 

/ [GV^ - ^VG] • dS = ^(ro, cj) + / GC'^dV , (9) 
Js Jv 

where S is the surface enclosing V, dS is directed out- 
wards and V = V — (A/'q /cj^) Vh • If 5* is the solar surface 
with radius i?, the expression is simplifled to 



or 



Jv 



(10) 

A flltered version of the left-hand side of Eq. (p!Q|) can 
be obtained from the observed wavefleld and a reference 
model. Moreover, if we compare the left-hand side of 
Eq. (fTOj) obtained from two sets of Green's functions, the 
differences will be due to the term with C . This includes 
what can be called the "anomalies" (e.g. magnetic ac- 
tivity), but also systematic errors in the wave operator 

As indicated in the introduction we will consider 
Green's functions G± corresponding to outgoing and in- 
going waves from the focus. In this case, Eq. (fTQ|) is a 
Kirchhoff representation for the wave fleld ^ in the fre- 
quency domain. We compute the egression and ingres- 
sion functions as 



^obs(^^,c^)^(^,c^;ro)d(^, (11) 
or 



where ^obs(^7 ^) is the observed signal (Doppler shift ve- 
locity in the frequency domain) and the integral is over a 
limited region of the solar surface, the pupil {V±). Cor- 
relating the ingression and egression functions one can 
expect to obtain a flltered version of the second term on 
the right-hand side of Eq. (fTQ|) . 

3. GRAVITO-ACOUSTIC RAY APPROXIMATION 

We will compare our results with those obtained from 
the ray theo ry. Rather than the aco ustic approxima- 
tion used by iLindsev fc BraunI (|2QQQb[ ) , we will consider 
a more general approximation. 

For high frequencies, the general solution of = 
can be expressed as a superposition of "rays" of the form 



^(r) = A{r) exp ( ±i jk dl 



(12) 



where we have introduced the local wave number, k. The 
integration is along a ray path and we have explicity 
indicated the two possible directions for each path. In 
this equation and hereafter we will removed the subindex 
for the reference model. 

In the ray theory, the ray path and the local wave vec- 
tor, k(r) = /c^r+kh, are determined from a flrst order 3D 
Liouville- Green expansion (a "WKB asymptotic approx- 
imation"). Speciflcally the dispersion relation is found 
to be given by 



... 0- (13) 

If the background state is spherically symmetric, k'^ = 
L'^/r'^, where L is a constant. 
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The ray path is determined from the group velocity. 
For the sphericahy symmetric case the rays are contained 
in a plane with a path given by 

^ = ^ . (14) 

dr rvr 

Here r and are the usual polar coordinates and the 
group velocity v = VrV -\-ve^ has components 

duj krUJ^c^ 



dkr - klc^N^ 



ve = 



khoo c 



(15) 



(16) 



The complex amplitude A(r) is obtained from a sec- 
ond order Liouville- Green expansion, called the trans- 
port equation. In our case this equation reduces to 



0. 



(17) 



Integrating this equation over the volume of a ray tube, 
applying the Gauss theorem and taking AS* for the 
tube section, one gets 



7V2 



kh ) AiS • n constant , 



(18) 



where n is a unit vector normal to the section. A rep- 
resentation of a ray tube si milar to that cons i dered h ere 
can be found in Figure 7 of lLindsev fc BraunI (|2QQQb[ ) . If 
we take the section such that n is in the radial direction, 
the amplitude can be determined by 



Ar(kr sin(9A(9) 



1/2 



constant . 



(19) 



where A^ is the angular size of the tube of section AS' 
0. 

3.1. Greenes functions 

The Green's functions defined by Eq. ([8]) have also 
a solution of the form ([12]) but with prescribed initial 
conditions, corresponding to a Dirac-delta impulse at the 
focus. Thus, the (complex) constant in Eq. ([19]) is fixed. 
Specifically, the amplitude for a Green's function at any 
point r far from the focus is given by 



^±(r|ro) 



exp (±z/co|r — re 



47r|re 



(kr sin6>)r 



kr sin 



1/2 



Ai9, 



1/2 



Ai9 
(20) 

Here re is a point close to ro where the integration along 
the ray path in ([T2]) starts up. For the Green's functions, 
here and in ([T2]) , the + sign corresponds to waves moving 
out of the focus ro (amplitude ^4+) and the — sign to 
waves moving into the focus (amplitude A_). 

We are interested in computing outward and inward 
Green's functions but, to be more general, after a given 
number of bounces at the surface and the interior. At the 
turning points, /c^ = and the amplitude cannot longer 
be computed by means of Eq. ffTO]) . In fact, we know 
that at the turning points ^(r) is just retarded by 7r/2. 
Since these Green's functions are defined only within the 
resonance cavity, both r and ro must be inside the cav- 
ity. In particular, at least formally, the Green's func- 
tions cannot be computed in the photosphere or above. 



Thus for a ray that initially goes inwards and then it is 
observed near the surface at its first arrival, the phase 
shift is — 7r/2. For any additional pair of inner and outer 
refiections the phase shift is — tt. So after s inner and 
s — 1 surface bounces (hereafter s skips), the phase shift 
is —(25 — l)7r/2 and the inward and outward Green's 
functions at a point r are given by 

.(2s- l)7r 



G±(r|ro) = A±(r|ro) exp > x 



expf±il f kdl^2{s-l) f kdl^ I kdlVk 

\ Wro^ri Jyi^Y2 Jyi^y ) / 



kdl\m) 



where ri and Y2 are the inner and upper turning points 
and yl±(r|ro) is given by (|2Q|) . The limits on the integrals 
indicate the interval over the ray path to be considered. 

4. SPHERICAL POLAR EXPANSION: ASYMPTOTIC 
APPROXIMATION 

For a spherically symmetric reference model, Eq. ([8]) 
can be solved by expanding the solution in terms of 
spherical harmonics in the form 



G(r|ro) = V/(r|ro) V Yl^{^o)YU^) = 



^=0 



47r ^ 

£=0 



g\r\ro){2i^l)P,{fi) (22) 



where Pei/J^) is a Legendre polynomial and fi = 
cos ^ cos ^0 + sin ^ sin 6^0 cos((/p — ^po). Here r ,6> and cp are 
standard spherical coordinates and Q = (^, cp). 

Since we are considering a spherically symmetric refer- 
ence model, we can choose, without losing generality, the 
polar axis in a such a way that ^o = and cpo = for the 
focus. In this way /i = cos 6 and the Green's functions 
have not dependence on (p. 

Inserting the decomposition given by (|22]) into Eq. ([8]), 
it is found that the radial factors geii satisfy the equation 



dr 



Kdi = at r / ro 



(23) 



where gi — rgi and kr is given by Eq. ([T3|) with = 
l{l + 1). The Dirac-delta term is transformed into the 
following inhomogeneous initial conditions: 



k{ro \ro) =gi{^o\^o) 



(24) 



and 



We also recall that the angular factors, given by the 
Legendre polynomials, satisfies the equation 



di92 



where Pi{cosO) = (csc^/^ l9)Q^((9) and 



^1 



1 1 



: sm 



(26) 



(27) 
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4.1. Asymptotic inwards and outwards Greenes functions 

We are interested in keeping the same kind of ingoing 
and outgoing Green's functions than in the ray theory, 
since this wih ahow us to do a straightforward compari- 
son. Since we have done a separation of radial and angu- 
lar variables, we need to search for inwards and outwards 
radial solutions and prograde and retrograde angular so- 
lutions and combine both in a suitable way. This kind 
of "travel- wave" solutions are well defined in the asymp- 
totic limit and we start with that. 

In a first order Liouville- Green expansion the general 
solution of the angular factors can be written as Q£{6) = 
(3e+Be^{0)^f3e-Be-{0), where 



Bi±{0) = ""^^^ exp [ ±i / n^dO 



(28) 



and PiJ^ and are constants. The value 6q will be 
taken at the focus. In a similar way the solution of the 
radial factors can be written as giir) = aiJ^CiJ^ir) + 
ai-Ci-{r) where and ai- are constants and 



Q±(r) 



-1/2 



exp 



±z / krdr 



(29) 



Here ri is the inner turning point. The constants and 
ai- can be obtained by introducing the equation gi{r) = 
a£-^Ci-^{r) -\- ai-C£-{r) into the initial conditions, 
and (|25]) . One gets 

<^£± = exp ( ^i J krdr^ . (30) 



2^0 y/kr{ro) 

In principle, we can take the outward prograde part 
of the solution, proportional to a^+(7^+(r)/3£+5^+(^), as 
the Green's function going out of the focus. However, 
the solution corresponds to an outgoing wave 

emanating from the focus ro and reaching the surface 
without any reflection. Since we are interested in the 
outgoing Green's function after s skips, we need to define 
new inwards and outwards asymptotic solutions. Let us 
call to the wave emanating inwards from the focus 
and now going outwards after s skips. This function and 
his inwards counterparts are: 



±i 



krdr - ^ 



(31) 
-7r/2 



where we have taken into account a phase jump of 
every time the wave reach a turning point. 

The angular factors has no turning points, so we 
can use the asymptotic functions defined in Eq. (|28)) . 
Then, the Green's function for a given degree / and 
frequency uj going out of the focus, proportional to 
ae+C!^{r)pe+Be+{e), is 

Gt+{r\ro,oj)-- 




The ingoing Green's function is the complex conjugate, 
e- - '^e+ ■ 



]e{0) = pe+{0)Be+{0) + (3,_{e)B,_{e) 



4.2. N on- asymptotic inwards and outwards Greenes 
functions 

Let us start with the angular factors. This corre- 
sponds to Eq. (|26|) that of course has the solution 
Qi = V^ndPi{cosO). To split it in prograde and ret- 
rograde parts we written the solution as 

(33) 

where B£± are still given by Eq. (|28|) . In the asymptotic 
approximation P£± were taken to be constants but now 
they are functions that change smoothly, at least where 
the asymptotic approximation is expected to work. Thus 
we replace the constants by osculating parameters that 
satisfies the additional condition 

Qm = (3e+{9)B',40) + /3e.{0)BUO) (34) 

where the primes means derivatives respect to 0. With 
this definition it is follow that 



{QiB[.-Q[Bij^) 



(35) 

Since Qi is known, j3i± and hence the prograde jS^^B^^ 
and retrograde parts Pi-Bi- are obtained. 

For reasons explained later, for the applications shown 
in this work we will use the asymptotic approximation 
for the radial factors. However, for completeness, we 
explain here how to obtain the non-asymptotic radial 
Green's functions. The general solution of Eq. (|23|) can 
be written as 



ge{r\ro) = a^+(r)C|+(r) + a^_(r)C|_(r) 



(36) 



where are given by (|3T]) . Here again, to fix the os- 
culating parameters we add the condition 

g',{r\ro) = a^+(r)(C|+)'(r) + a,_(r)(C|_)'(r) (37) 

from which it is found that 



ae+ = -mCW-g'eCl 



at- = --{m{Ch)'-~9[Ct+). 

(38) 



Thus, for a given solution ^£(r|ro), the functions a^+(r) 
and a^-ir) can be computed. Note that gi{r\rQ)] satisfies 
Eq. (|23|) and since we are dealing with waves reflected 
in the interior, the solution must be regular at the cen- 
ter. With this condition alone, ^£(r|ro) is determined up 
to a constant factor. On the other hand ^^(ro|ro) can 
be taken as the one given by the asymptotic solution, 
namely ^^(ro) = Q^^-(ro)C|_ (ro) + a£+(ro)C|^(ro) , with 
Oii±{rQ) given by Eq.([^ and C^i^ by Eq.([3lD. This fixes 
^^(^ko) completely and hence the inward and outwards 
solutions, a£-{r)Cf_{r) and a£+(r)C|^(r) respectively. 
Combining the outward solution with the prograde so- 
lution of the angular part, the outgoing wave is given 

by 



G|+(r1fo) 



2? + 1 
167rr 



csc^/2^ X 



{Qi{e)B[_{9) - Q'i{e)Bi_{e)} b,+(^) x 

{Ur)CUr) - g[{r)C,-{r)] C,+ (r) . (39) 
The ingoing Green's function is the complex conjugate. 
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Fig. 1. — Phase shift computed using equation (|42|) . Eigen- 
frequencies, with degrees i = 20 — 80, are taken from MDI data. 
Model S of Christensen-Dalsgaard (1996) have been used. The 
sohd thick hue represent a fitted polynomial of degree 10. 

5. THE PHASE SHIFT 

In the uppermost layers of the Sun, the asymptotic ap- 
proximation fails. Moreover, non adiabatic effects and 
the dynamics of convection, that are not included in 
the wave equations, are too large in these layers to be 
neglected. These effect s were taking into account by 
iLindsev fc BraunI (|2QQ4D by using empirical control cor- 
relation functions. Here we use a completely different 
approach. 

When computing the wave functions in the asymptotic 
approximation, these effects can be taken into account 
by introducing an additional phase shift at t he upper 
turning point — e.g. Gough (1984), iDeubner fc Gough 
(|1984) — . In the same way the Green's functions can be 
modified to include such effects, namely 

gi{r,0\ro,Oo) = Gt{r,0\ro,Oo) exp {2isn5a) (40) 

where 27rda — 7r/2 is the phase shift after one surface 
reflection. To a good approximation, Sa is a function 
of frequency alone. If we were to include just the formal 
bounces at the surface, s = s — 1. However the pupil and, 
depending on the problem, the focus are so close to the 
upper turning points that each of them can account for 
an extra shift in the phase of, say, a half bounce. Thus 
in practice one should take s = s — 1/2 if the focus is 
below the uppermost layers or s = 5 if the focus is too 
close to the surface^ 

As shown by Ch ristensen-Dalsgaard fc Perez HernandezI 
(|T992l ). a theoretical phase shift at every point, 6a{r)^ 
can be computed by comparing the asymptotic ap- 
proximation and the numerical solutions of the wave 
equations. This would allow to use Eq. (|4Q|) to correct 
the Green's functions at any point in the interior, but, 
unfortunately, our tests have shown that this correction 
is not enough to achieve good seismic images. The 
reason is that they did not include any non adiabatic or 
convective terms in the full wave equations, and this is 
not a simple task. 

For this work, we have computed an r-independent 
phase shift based on the observed frequencies of the nor- 
mal mode of oscillations. This way we include all the 
neglected aspects in the wave equations but the focus 
must be assumed to be above or below these uppermost 
layers. Moreover, the method is limited to the asymp- 
totic approximation for the radial factors of the solution. 
Specifically, 5a{uj) can be estimated by inserting the ob- 



served frequencies into the asymptotic expression 

J krdr = TT - ^ - Sa^ , (41) 

where n is an integer, the radial order of the mode. In 
our case the wave number is given by Eq. (p!3|) and hence 




(42) 

To obtain fo(a;), eigenfrequencies from u = 2500 /iHz 
to = 4500 /iHz and degrees = 20 — 80 were taken 
from MDI data. Modes of lower degrees are rejected be- 
cause they do not satisfies the Cowling approximation as 
good as intermedia te and high-degree modes. In addi- 
tion. Model S from IChristensen-Dalsgaard et al.l (| 19961 ) 
was used. A fitted polynomial of degree 10 to Eq. (j42j) 
was considered. The result is presented in Fig. [H Notice 
that TvSa changes by more than 1.5 r ad in the frequency 
range of interest. 

6. NUMERICAL RESULTS 

To validate the Green's functions calculated using this 
new approach, we apply them to create seismic maps of 
active regions at the far side of the Sun. In the case of 
far-side seismic holography, the focus is located at (or 
just below) the surface of the far hemisphere of the Sun 
and the pupil is an annulus surrounding the antipodes 
of the focus. Although seismic maps of the full far-side 
hemisphere can be calculated using waves following one, 
two and three-skip ray paths before arriving at the pupil 
(i Brau n fc Lindsey 2001), for simplicity, here we will con- 
centrate only on the central part of the maps. Tradition- 
ally, this is calculated consid ering waves following a 2 - 
skip ray path, as described bv lLindsev fc BraunI (|2000b[ ). 

6.1. Comparison of ray theory and polar expansion 

As it has been done in previous work, we normalize the 
Green's functions at each frequency to its mean absolute 
value. Fig. [2] shows these normalized Green's functions 
against the polar angle for different frequencies. Results 
for the polar expansion and the gravito-acoustic rays are 
compared. Eq. (|39|) has been used for the polar ex- 
pansion, but with the asymptotic approximation for the 
radial factors, while for the ray theory Eq. (|2T]) has been 
considered. In both cases the surface phase shift shown 
in Fig. [T]has been added according to Eq. (|4Q|) . Finally, 
what we actually show are the derivatives of the Green's 
functions that we will use to compute the ingression and 
egression fields following Eq. (pTj) . The focus is located 
at ro = O.9995i?0, = and the Green functions are 
shown at the upper turning point after two skips. 

In the ray theory, for a given number of skips, and a 
given frequency cj, each point in the pupil corresponds 
to a given value of the horizontal wave number at the 
focus, kh{ro) and hence of L = i{i + 1). Here ^ is a 
real constant. Our pupil goes from 6>min = 2 to ^max = 
3 rad; this corresponds to and interval in "degrees" from 

^min — 15 to /max — 28 for U = 2500 /iHz and Imin — 27 to 

/max — 50 for u = 4500 /iHz. In the polar expansion there 
is not such relation, in fact in Fig. [2] Green's functions 
for degrees up to / = 60 have been added up, the same 
for all the frequencies. 
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Fig. 2. — Comparison of Green's functions at selected frequencies. 
The black continuous lines correspond to the Green's functions 
calculated using the Legendre polynomial decomposition and the 
red dashed ones to those calculated using a gravito-acoustic ray 
path approximation. 

As it can be seen in Fig. [21 there is a systematic shift 
between the Green's functions computed with the two 
different approaches, but the solutions are closer at low 
angles; in fact in the ray theory these angles correspond 
to higher / where this approximation is expected to work 
better. On the other hand the solutions of the polar 
expansion approximation look more irregular at low fre- 
quencies; this is related to the fact that at these frequen- 
cies the upper turning point r2 is too close to the focus 

We use both sets of Green's functions described in this 
paper, as well as the traditional ones (acoustic ray ap- 
proximation), to calculate maps of magnetic activity at 
the non- visible hemisphere of the Sun. Fig. [3] presents 
far-side maps calculat ed using; Globa l Oscil lation Net- 
work Group (GONG, iHarvev et aTl data for 
September 1 and 2 2005 for the three sets of Green's 
functions. The maps clearly show the signature of active 
region NOAA- 10808 as a large negative phase shift. The 
Green's functions used for the analysis include frequen- 
cies computed from u = 2500 /iHz up to 4500 /iHz with 
a step of ~ 8 /iHz. Angles go from 6 = 2 to 3rad with a 
step of 0.0016 rad. Active region NOAA 10808 appeared 
at the East limb of the Sun September 7 2005. 

The signature of the active region is very similar in 
the maps for all the cases. However, while the noise level 
is similar in all three cases, the spatial distribution is 
only correlated between the maps calculated using the 
two types of Green's functions described in this paper. 
This may be to the different approaches used for the 
sub-surface phase shift correction, explained in section 
5. However more research is necessary in order to under- 
stand the different contributors to the noise. 

6.2. Green's functions combining different number of 

bounces 

As mentioned before, in the polar expansion approx- 
imation there is not a one to one relation between the 
angle in the pupil and the mode degree. Thus, even for 
the same pupil and the same angular degrees it is possible 
to combine the Green's functions from different bounces. 
In Fig. Hlwe show the Green's functions from the polar 
expansion against the polar angle. In this case we have 
used degrees from ^ = 15 to ^ = 100. The Green's func- 
tions are shown for 5 = 2, s = 3 and the sum of both 
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Fig. 3. — Post el-projected maps of the far side calculated us- 
ing the traditional acoustic ray approximation (top), the gravito- 
acoustic ray approximation (center) and the spherical polar ex- 
pansion (bottom) Green's functions. The maps extend to approxi- 
mately 60 degrees from the antipodes of the center of the solar disk 
facing the Earth. The areas of large negative phase shift (cj)) are 
related to areas of strong magnetic field. In this case, the strong 
signature corresponds to active region NOAA 10808 5-6 days before 
it appeared at the East limb of the Sun on September 7 2005. 
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Fig. 4. — Normalized Green's functions against polar angle from 
the focus computed with the polar expansion technique for s = 2, 
s = 3 and the sum of both. Black lines are for v = 3000 /i-,Hz and 
the red ones for v = 4500 /i,Hz. 

(with the normahzation at each frequency done after the 
sum). We do not consider the s = 1 case since most of 
the signal for the first bounce in our pupil comes from 
very low degrees for which the Cowling approximation 
does not work properly. As it can be seen, the maximum 
amplitude of these Green functions change both with s 
and frequency. 

We calculate a series of far-side maps using the polar- 
expansion Green's functions that combine two and three 
skips. In Fig. [5l we compare these maps with those using 
Green's functions that include only 5 = 2. A sequence 
of three days spanning November 27 to 29 2006 show- 
ing a weaker active region just south of the equator an 



7 



Nov. 27 2006 Nov. 28 2006 Nov. 29 2006 




X (rodions) X (rodions) X (rodions) 



Fig. 5. — Postel projected maps of the far side calculated using 
the spherical polar expansion Green's functions for two skips (top) 
and a combination of two and three skips (center). The signature 
of active regions NOAA 10930 can be seen as at the far side hemi- 
sphere, it appeared at the front side on December 5 2006 just below 
the solar equator. The bottom panels show a cut of the far-side 
maps at heliocentric latitude -2.3 degrees. The phase shift associ- 
ated to the magnetic area is enhanced when using the combination 
of 2 and 3 skips by approximately 20-30% compared to the 2 skips 
approach used for the center of the map. 

past the central meridian is presented. This active region 
was at the same location in the front site the previous 
Carrington rotation (NOAA 10923) and reappeared at 
the East limb December 5 2006 (NOAA 10930). We find 
that the signature associated to the magnetic area is en- 
hanced when using the combination of 2 and 3 skips by 
approximately 20-30% compared to the 2 skips approach 
used for the center of the map. 

7. CONCLUSIONS 

We have implemented a holographic, or migration tech- 
nique, that includes a factorization of the Green's func- 
tions in angular and radial parts. This technique is more 
general than the one provided by the ray theory. We 
show examples of how these new Green's functions can 
be used to improve far-side seismic images by combining 
the Green's functions corresponding to travel waves with 
t wo and three skips or inner reflec t ions. 

IZh^ ([2007h and lllonidis et all using a time- 

distance technique have demonstrated that adding more 
skips to the calculation of far-side maps improves the 
signal-to- noise. However, in their analysis, they combine 
seismic images obtained from different skip ray paths af- 
ter they are calculated. Here, we show that by combin- 
ing the Green's functions associated to two and three 
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skips before calculating the map in seismic holography 
we improve the signal to noise in a similar way. Since 
the spherical geometry of the observational pupil makes 
the calculation of the ingression and egression functions 
for far-side maps computationally expensive, the solution 
presented here is an ideal way of improving the signal to 
noise without the need for increasing the computation 
time. Our approach is also formally different, since we 
sum the Green's functions from different bounces with- 
out a previous normalization. 

There are some previous works where Green's functions 
beyond the ray approxim ation were used, for instance in 
iLindsey fc BraunI (|2004l ) , but there a completely piano- 
parallel approximation was considered and hence the 
technique was limited to focuses in the front-side and 
close to the surface. On the other hand for our Green's 
functions, although we have shown here the application 
to the particular case of far-side seismic holography, the 
technique is general enough to be applied in the front side 
and for focuses under the surface, even below the convec- 
tion zone, at least formally. In that scenario, the signal- 
to-noise increases considerably and more work would be 
required to test the inferences. 

An interesting result of our work is the fact that the 
sub-surface phase shift, that has been recognized as a 
fundamental correction to any Green's function in the 
past, can be computed by using the observed frequencies 
of the normal modes of oscillations. This opens a link be- 
tween local and global helioseismology and perhaps the 
issues of the local approach can be used to learn some- 
thing about the global one and vice versa. 
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